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Bayesian statistical inference with a case-crossover design was used to examine the effects of air pollutants 
{Particulate matter <10 urn in aerodynamic diameter (PM 10 ), sulphur dioxide (S0 2 ), and ozone (0 3 )} on 
mortality. We found that all air pollutants had significant short-term impacts on non-accidental mortality. 
An increase of 10 ug/m 3 in PM 10) 10 ppb in 0 3 , 1 ppb in S0 2 were associated with a 0.40% (95% posterior 
interval (PI): 0.22, 0.59%), 0.78% (95% PI: 0.20, 1.35%) and 0.34% (95% PI: 0.17, 0.50%) increase of 
non-accidental mortality, respectively. 0 3 air pollution is significantly associated with cardiovascular 
mortality, while PM 10 is significantly related to respiratory mortality. In general, the effects of all pollutants 
on all mortality types were higher in summer and winter than those in the rainy season. This study highlights 
the effects of exposure to air pollution on mortality risks in Thailand. Our findings support the Thailand 
government in aiming to reduce high levels of air pollution. 



As Thailand transfers its economic growth from an agricultural to an industrial base 1 , it faces increasing 
levels of air pollution which have been found in other countries to be associated with adverse health 
consequences. Particulate matter <10 urn in aerodynamic diameter (PM 10 ) is the most important air 
pollutant in both urban and rural areas, the concentration of which is as high as or higher than those in North 
America and European countries 2 . Meanwhile, Thailand is also experiencing high ozone (0 3 ) pollution because 
of the intense emission of air pollution coupled with year-round sunlight 3 . Even in the rural areas, forest fires and 
agricultural burning have contributed to high levels of air pollution, which have increased to critical levels since 
2006 1,4 . However, in order to continue the economic growth, expansion of petrochemical plants and industries 
rose sharply, which greatly contribute to air pollution. 

Many studies worldwide have shown that short-term exposure to air pollution is associated with mortality and 
morbidity including cardiovascular and respiratory diseases, exacerbation of chronic respiratory conditions, and 
decreased lung function 5 9 . However, only a few studies on air pollution and mortality have been conducted in 
Thailand 210 . Data were limited to Bangkok, where the estimated effect of air pollution on mortality was higher 
than Hong Kong, Shanghai and Wuhan 10 . Measurements of air pollution-health associations are limited in other 
areas in Thailand. In addition, single-city studies have limited generalisability and it is difficult to combine their 
estimates of air pollution effects, because different modelling approaches have been used. Therefore, it is neces- 
sary and important to conduct a multi-city analysis of air pollution in Thailand, where the characteristics of air 
pollution (e.g., level and components) and socio-demographic status of local residents (e.g., disease pattern, and 
socioeconomic status) are different from those of developed countries. 

In this study, we examined the association between daily air pollutants (PM 10 , 0 3 , and sulphur dioxide (S0 2 )) 
and mortality (non-accidental, cardiovascular and respiratory diseases) in 18 Thailand provinces (Figure 1). We 
also examined whether the associations differed by season. 

Results 

Table 1 summarises mortality, air pollution, temperature and relative humidity data in the 18 Thailand provinces. 
Daily average counts of non- accidental, cardiovascular, and respiratory deaths differed by province and ranged 
from 4 to 66, from 1 to 14, and from 1 to 8, respectively. Daily mean concentrations of PM 10 , 0 3 and S0 2 in the 
Thailand provinces ranged from 31.2 to 94.9 ug/m 3 , from 9.9 to 25.6 ppb, and from 1.1 to 13.9 ppb, respectively. 
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Figure 1 | The mean values of (A) non-accidental death, (B) PM 10 , (C) 0 3 , and (D) S0 2 in 18 Thailand provinces during 1999-2008. R software with 
"maptools package" was used to produce the maps. 



All provinces had similar temperature and relative humidity. In gen- 
eral, the air pollution levels are higher in north middle provinces than 
southern Thailand (Figure 1). Similar pattern was found for non- 
accidental deaths. 

Table 2 shows the mean values for the Spearman correlation coef- 
ficients between air pollutants and weather conditions in 18 pro- 
vinces in Thailand. PM 10 , 0 3 and S0 2 were positively correlated 
with each other, except that 0 3 and S0 2 had slightly positive cor- 
relation. Temperature and relative humidity were negatively corre- 
lated with PM 10 and 0 3 . The results for each province are shown in 
supplemental material Table SI. 

Figure 2 shows the pooled percentage increase (mean and 95% 
Pis) of mortality associated with air pollutants, using single-pollutant 
models with different lag structures in Thailand during 1999-2008. 
All air pollutants had short-term effects on mortality. In general, for 
all death categories, the effect estimates achieved highest at lag 0 for 
PM 10 , at lag 0-5 for 0 3 , and at lag 0-3 for S0 2 . Thus, we used these 
lags for each pollutant in the following analyses. 



Figure 3 shows the percentage increase (mean and 95% Pis) of 
mortality associated with air pollutants in 18 Thailand provinces 
during 1999-2008. It is clear that the effect estimates of air pollutants 
vary greatly by provinces. We observed statistically significant 
impacts of PM 10 , 0 3 and S0 2 on non-accidental, cardiovascular, 
and respiratory mortality in most provinces. Heterogeneity for all 
associations between air pollutants and mortality was significant 
(results not shown). 

Table 3 shows the pooled percentage increase (mean and 95% Pis) 
of mortality associated with air pollutants in Thailand during 1999- 
2008, using single and multi-pollutant models. In the single pollut- 
ant models, an increase of 10 ug/m 3 in PM 10 , 10 ppb in 0 3 , 1 ppb in 
S0 2 were associated with a 0.40% (95% PI: 0.22, 0.59%), 0.78% (95% 
PI: 0.20, 1.35%) and 0.34% (95% PI: 0.17, 0.50%) increase of non- 
accidental mortality, respectively. In the two/three-pollutant models, 
the estimated associations between air pollutants and mortality were 
attenuated in comparison with single-pollutant models, but gen- 
erally the associations between air pollutants and non-accidental 
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Table 1 | Descriptive summary for the study period, daily deaths, air pollutants and weather conditions in 1 8 Thailand provinces, during 
1999-2008 

Mean (Standard deviation) 



Non-accidental Cardiovascular Respiratory PMio Mean Relative 

Province Period death death death (ng/m 3 ) O3 (ppb) S0 2 (ppb) temperature (°C) Humidity (%) 
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mortality remained significant. 0 3 is significantly associated with 
cardiovascular mortality, while PM 10 is significantly related to res- 
piratory mortality. 

Table 4 presents the pooled effect estimates (mean and 95% Pis) 
for the increase in mortality associated with air pollutants in three 
seasons in Thailand during 1999-2008. The associations between air 
pollutants and mortality differed by season. We observed that the 
associations between air pollutants and mortality were higher in 
summer and winter than those in rainy season. For example, an 
increase of 10 |ig/m 3 in PM 10 was associated with a 0.51% (95% PI: 
0.28, 0.74%), 0.20% (95% PI: -0.09, 0.48%) and 0.35% (95% PI: 0.05, 
0.65%) increase of non- accidental mortality in summer, rainy season 
and winter, respectively. 

Sensitivity analyses show that the change of time strata and intro- 
ducing N0 2 and CO did not substantially affect the associations of air 
pollutants with mortality (results not shown), suggesting that our 
findings are relatively robust in this aspect. However, when we 
removed temperature and relative humidity from the models, the 
effect estimates of air pollutants on mortality were higher than those 
without controlling for temperature and relative humidity 
(Supplemental Material, Table S2). This suggests that it is necessary 
to control for temperature and relative humidity when assessing the 
association between air pollution and mortality. When we used 
empirical statistical models (generalised linear model and general- 
ised additive model) for this study, the mean effect estimates are 
similar as those from the Bayesian statistical models. However, the 
confidence intervals from empirical models are wider than those 
from Bayesian models. 



Table 2 | Average values for Spearman correlation between air 
pollutants and weather conditions in 1 8 Thailand provinces 



Pollutant O3 SO2 Mean temperature Relative humidity 



PM 10 0.45 0.18 -0.17 -0.39 

O3 0.03 -0.23 -0.52 

SO z 0.01 -0.04 

Mean temperature 0.03 



Discussion 

We used Bayesian statistical inference with the case-crossover design 
to examine the associations between air pollutants and mortality in 
18 provinces in Thailand. To the best of our knowledge, it is the 
largest study to systematically assess the effects of air pollution on 
mortality in Thailand. We found that all pollutants (PM 10 , 0 3 , and 
S0 2 ) had significant impacts on non-accidental mortality when we 
used single pollutant models. But the magnitude of effect estimates of 
air pollutants on mortality were attenuated when two- or three-pol- 
lutant models were used. The effects of air pollutants on mortality 
appeared immediately and lased for 2 days for PM 10 , 5 days for 0 3 , 
and 3 days for S0 2 . The effects of air pollutants on mortality were 
higher in summer and winter than those in rainy season. 

For PM 10 . In general, our effect estimates associated with PM 10 in 
Thailand are consistent with previous studies, except for 
cardiovascular mortality. For example, our estimate for non- 
accidental mortality is 0.40% (95% PI, 0.22-0.59%). A study using 
75 single-city data worldwide produced an effect estimate of 0.6% 
(95% confidence interval (CI), 0.5-0.7%)". An USA study with 90 
cities generated an estimate of 0.2% (95% CI, 0.1-0.4%) 12 , while an 
European study of 29 cities gave an estimate of 0.6% (95% CI, 0.4- 
0.7%) 13 . A study in 14 USA cities generated an estimate of 0.35% 
(95% CI, 0.2-0.5%) 14 . A multi-city study in four Asian cities 
(Bangkok, Seoul, Inchon, and Hong Kong) gave a pooled estimate 
of 0.55% (95% CI, 0.26-0.85%) 10 . 

However, some higher estimates have been reported in other 
developing countries, for example, a Mexico City study showed that 
an effect estimate of 1.8% (95% CI, 0.9-2.7%), while a result for 
Santiago, Chile was 1.1% (95% CI, 0.9-1.4%) 15 ' 16 . A previous study 
conducted in Bangkok using data from 1999-2003 produced an 
estimate of 1.3% (95% CI, 0.8 to 1.7%), which was higher than our 
result in Bangkok. The reasons might be that previous studies in 
Bangkok and other developing counties used relatively old data. In 
addition, it may be that people have recently changed their behaviour 
to protect themselves from serious air pollution, for example, using 
mouth mask during the heavily polluted days. 

In this study, we found that the estimated effects of PM 10 on non- 
accidental mortality were similar to respiratory mortality and were 
higher than for cardiovascular mortality. This is contrary to our 
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Figure 2 | Percentage change (mean and 95% Pis) of mortality associated with a 10 (ig/m 3 increase in PM^ 
in S0 2 , using different lag structures in Thailand during 1999-2008. 
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Figure 3 | Percentage change (mean and 95% Pis) of mortality associated with a 10 ng/m 3 increase in PM 10 (lag 0), 10 ppb increase in 0 3 (lag 0-5), and 
1 ppb increase in S0 2 (lag 0-3) in 18 Thailand provinces during 1999-2008. 
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Table 3 | Pooled effect estimates (mean and 95% Pis) for the increase in mortality associated with an increase of 1 0 ug/m 3 in PMio, 1 0 ppb 
in O3, 1 ppb in SO2 in Thailand during 1 999-2008, using single-/multi-pollutant models 

% increase (Pis) 

Model Non-accidental Cardiovascular Respiratory 



PM 10 

Adjust for O3 


0.40 (0.22, 0.59) 


0.13 (-0.16, 0.41) 


0.40 (0.07, 0.73) 


0.37(0.18, 0.56) 


0.09 (-0.18, 0.37) 


0.32 (-0.02, 0.66) 


Adjust for S0 2 


0.37(0.20, 0.55) 


0.09 (-0.19, 0.38) 


0.34 (0.01,0.68) 


Adjust for O3 and SO2 


0.35 (0.17, 0.52) 


0.07 (-0.21,0.34) 


0.28 (-0.06, 0.62) 


O3 


0.78 (0.20, 1.35) 


1.25 (0.15, 2.36) 


0.93 (-0.65, 2.51) 


Adjust for PM10 


0.49 (-0.05, 1.03) 


1.20 (0.16, 2.25) 


0.66 (-0.72, 2.04) 


Adjust for S0 2 


0.64 (0.05, 1.23) 


1.12 (-0.03, 2.27) 


0.69 (-0.85, 2.23) 


Adjust for PM10 and S0 2 


0.40 (-0.1 6, 0.96) 


1.09 (-0.01, 2.19) 


0.47 (-0.91, 1.86) 


so 2 


0.34 (0.17, 0.50) 


0.33 (-0.04, 0.70) 


0.49 (-0.01, 0.98) 


Adjust for PM10 


0.21 (0.08,0.41) 


0.30 (-0.07, 0.67) 


0.38 (-0.09, 0.86) 


Adjust for O3 


0.31 (0.14,0.47) 


0.27 (-0.10, 0.64) 


0.44 (-0.05, 0.93) 


Adjust for PM10 and O3 


0.23 (0.06, 0.40) 


0.26 (-0.12,0.64) 


0.36 (-0.1 1, 0.83) 



hypothesis that those with cardiopulmonary disease should be more 
sensitive to air pollution than non-accidental mortality. Many pre- 
vious studies have confirmed that mortality risks from PM 10 were 
higher for cardiovascular and respiratory mortality 1017 . The reason 
might be that the effects of PM 10 on mortality varied greatly by 
province (Figure 3), which makes the pooled effect non-significant 
for cardiovascular mortality. In addition, the low effect values esti- 
mated by our model could be due to publication bias in other single- 
city studies that are incorporated into the meta-analyses. Because we 
used the same statistical approach to the 18 large Thailand provinces, 
our results are free from publication bias. However, the results in 
Bangkok were similar to previous studies which reported that PM 10 
had significant impacts on cardiovascular and respiratory mortal- 
ity 2,10 . Further detailed studies are needed to explore the reason. 

For 0 3 . Our results indicate that 0 3 has a substantial health burden 
on mortality in Thailand. According to our pooled national estimate 
from the single pollutant model, a 10 ppb increase in daily 0 3 was 
associated with an increase of 0.78% (95% PI: 0.20, 1.35%) in non- 
accidental mortality. The effect estimates for cardiovascular and 
respiratory mortality were slightly higher than the one for non- 
accidental mortality. Our results are consistent with previous 
studies 18 . Meta-analysis study showed that the pooled estimates 
expressed as the percentage increase in mortality for a 10 ppb 
increase in daily 0 3 were 0.87% (95% CI, 0.55%-1.18%) for non- 
accidental mortality, 1.11% (95% CI, 0.68%- 1.53%) for 
cardiovascular mortality; and 0.47% (95% CI, -0.51%-1.47%) for 
respiratory mortality 19 . 

Ozone pollution is still a problem in both urban and rural areas in 
Thailand 20 , and is related to health problems as we found in this 
study. In urban areas, its rise was caused primarily by increased usage 



of motor vehicles, wherein vehicle emissions are a major source of 
precursor hydrocarbons and nitrogen oxides. In rural areas, agricul- 
tural burning and forest fires have contributed to high levels of air 
pollution, which have increased to critical levels since 2006 1,4 . This 
increased biomass burning could potentially enhance tropical tro- 
pospheric ozone. At the same time, the higher UV intensities com- 
bined with higher humidity in Thailand's tropical atmosphere result 
in an increased amount of O3 3 . Our results, combined with previous 
epidemiologic and toxicological evidence on ozone toxicity, suggest 
that this widespread 0 3 pollution may adversely impact human 
health. 

For S0 2 . This study showed evidence that outdoor S0 2 was related to 
the increased risk of non-accidental mortality in Thailand. The effect 
estimates of S0 2 on cardiovascular and respiratory mortality were 
positive but non-significant. Our findings were generally comparable 
with most previous studies. Our results in Bangkok were similar to a 
previous study 21 . A meta-analysis study reported that the excess non- 
accidental mortality risk associated with a 10 ug/m 3 increase of S0 2 
was 1.00% (95% CI, 0.75 to 1.24%) in four Asian cities (Bangkok, 
Hongkong, Shanghai and Wuhan) 21 . Another meta-analysis 
indicated that a 10 ug/m 3 increase of S0 2 was associated with a 
0.75% (95% PI, 0.47 to 1.02) increase of non-accidental mortality 
in 17 Chinese cities 22 . However, some North American studies did 
not report a significant association between S0 2 and mortality, which 
is possibly because of the very low levels of S0 2 air pollution 23,24 . 

The estimated mortality risks associated with S0 2 in Thailand 
provinces were consistent with previous studies worldwide. 
However, the burden of disease from exposure to S0 2 might be lower 
than that in other countries because the S0 2 levels in Thailand are 
much lower than those reported in other developing countries for 



Table 4 | Pooled effect estimates (mean and 95% Pis) for the increase in mortality associated with an increase of 10 ug/m 3 in PM 10 , 1 0 ppb 
in O3, 1 ppb in SO2 in three seasons in 1 8 Thailand provinces during 1 999-2008, with single-pollutant model 

% increase (Pis) 

Pollutant Season Non-accidental Cardiovascular Respiratory 



PM 10 


Summer 


0.51 (0.28, 0.74) 


0.54 (0.03, 1.05) 


0.58 


-0.05, 1.20) 




Rainy 


0.20 (-0.09, 0.48) 


-0.51 (-1.4,0.38) 


0.24 


-0.56, 1 .04) 




Winter 


0.35 (0.05, 0.65) 


-0.02 (-0.56, 0.51) 


0.33 


-0.12, 0.78) 


O3 


Summer 


0.52 (-0.1 5, 1.20) 


1.92 (0.39, 3.45) 


1.52 


-0.34, 3.38) 




Rainy 


0.37 (-0.73, 1.46) 


0.37 (-2. 10, 2.84) 


1.21 


-2.98,5.40) 




Winter 


1.23 (0.27,2.19) 


1.56 (-0.1 1, 3.24) 


0.79 


-1.96, 3.55) 


so 2 


Summer 


0.32 (0.02, 0.62) 


0.71 (0.007, 1 .40) 


0.48 


-0.33, 1.30) 




Rainy 


0.08 (-0.25,0.41) 


0.15 (-0.60, 0.91) 


0.22 


-0.94, 1.38) 




Winter 


0.57 (0.31,0.84) 


0.31 (-0.048, 0.87) 


0.72 


-0.14, 1.58) 
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example, China 21,22 . Also, the population in Thailand is lower than in 
China. However, ambient S0 2 may still represent a major public 
health concern in Thailand, as we found it has a significant effect 
on non- accidental mortality. These findings contribute to the scient- 
ific literature on health effects of S0 2 for low exposure settings typical 
in developing countries. 

Lagged effects. Understanding the short-term lag pattern of the 
association between air pollution exposure and mortality is very 
important for healthcare providers and public health authorities in 
development of response plans. This study showed that the impacts 
of air pollution on mortality risks were mainly limited to recent days 
of exposure (lag 0-5 days) in Thailand. Our findings are consistent 
with previous studies of air pollution-mortality relationships 22 ' 25 , 
which suggest that the effects of air pollution peak quickly. 
Therefore, applying timely preventive measures are effective in 
reducing the health effects of air pollution in Thailand, when 
extreme air pollution days appear. However, this will depend on 
timely responses by the health alert system, and better awareness 
of the dangers of air pollution by government, public health 
professionals, medical professionals and the general public. 

Multi-pollutant models. We found all pollutants (PM 10 , 0 3 and 
S0 2 ) have a great impact on mortality risks (except PM 10 and 
cardiovascular mortality) in Thailand, when we used single 
pollutant models. Therefore, strict standards for air pollutants 
should be considered when establishing Thailand's air-quality 
strategies. When we used the multi- pollutant models, the 
magnitude of effect estimates of air pollutants were attenuated. 
These findings were consistent with previous studies 17,22 ' 25 . This 
might be caused by the co-linearity between air pollutants, but it is 
unlikely from the viewpoint of biological mechanisms. However, it is 
very difficult to separate the effects in multiple pollutants models. 
Our results of two or three pollutant models should be interpreted 
with caution, because both S0 2 and 0 3 are precursors of secondary 
particles. Thus, our ability to determine precisely the independent 
effects of individual pollutants on the risk of death is limited. 

Seasonal modification. We found strong evidence that PM 10 , 0 3 
and S0 2 had significant effects on mortality in Thailand. In 
general, these estimated effects were higher in summer and winter 
than those in rainy season. Similar findings have been reported by 
previous studies 26 ' 27 . The main reason is that in the rainy season 
frequent rainfall reduces air pollution (Supplemental Material, 
Table S3-S5). In addition, in the winter and summer, the effects of 
air pollution on mortality might be enhanced by both cold and hot 
temperatures 28,29 . 

Our study has some limitations. We used ambient pollutant con- 
centrations as individual exposure, which might induce unavoidable 
measurement error. The local air monitoring stations are located in 
urban areas at least 50 meters distant from a main road. This may 
introduce measurement bias to rural areas. Because the monitoring 
stations are located beside the main roads, it is potential to introduce 
biases, for example, upwind/downwind of major point sources or 
highways would influence the individual exposure. The data are at 
province level, so the effects of air pollutants on mortality in Thailand 
may be underestimated. 

In conclusion, there is clear evidence that air pollution can affect 
existing health conditions and greatly increase mortality 10 . In addi- 
tion, air pollution shortens life expectancy and health expectancy 30,31 . 
Thus, it is necessary to pay attention to the dangers of air pollution 
and to inform the public on how to minimise their risks from serious 
air pollution. This study highlights the effects of exposure to air 
pollution on mortality risks in Thailand. Our findings support the 
Thai government efforts in reducing high levels of air pollution in 
Thailand, in order to protect the health of the population. 



Methods 

Data collection. This study includes 18 Thailand provinces: Ayutthaya, Bangkok, 
Chachoengsao, Chiang mai, Chon buri, Khon kaen, Lampang, Nakhon sawan, 
Nakhon ratchasima, Nonthaburi, Pathum thani, Ratchaburi, Rayong, Samutprakan, 
Samut sakhon, Saraburi, Songkhla, Surat thani (Figure 1). Our study areas include 
both urban and sub-urban areas of these provinces, because both are experiencing the 
same air pollution issues in Thailand and Thai provinces are similar in area to big 
cities in the USA or China. In addition, death counts, and therefore statistical power, 
are limited if we restrict the study to urban areas. Previous studies have used similar 
data to examine the effects of air pollution on mortality in China 32 . 

We obtained daily counts of non- accidental death in 18 Thailand provinces from 
the Ministry of Public Health, Thailand during 01/011999 to 31/12/2008. Causes of 
death were coded according to the International Classification of Diseases, Tenth 
Revision (ICD-10). We classified the mortality data as deaths due to non-accidental 
causes (ICD-10: A00-R99), cardiovascular disease (ICD-10: 100-199), and respir- 
atory disease (ICD-10: J00-J99). 

We obtained daily data on air pollution from the Pollution Control Department, 
Ministry of Natural Resources and Environment for 18 provinces during 1999-2008 
(Air Quality and Noise Management Bureau 2010). Monitoring air pollution data 
started from 2004 for Chachoengsao, and from 2000 for Nakhon sawan. For each 
province and air pollutant, daily concentration was averaged by fixed air quality 
monitoring stations within the province. If monitored data for an individual pollutant 
were insufficient to calculate a daily average, all measurements from that day were 
excluded for that pollutant and monitor. The averages of daily mean concentrations 
for five air pollutants were obtained in this study: 24-h average of PM 10 , 24-h average 
of ozone (0 3 ), 24-h average of S0 2 , 24-h average of nitrogen dioxide (N0 2 ), and 8-h 
average of carbon monoxide (CO). We only examined the effects of PM 10 , 0 3 and S0 2 
on mortality in this study, because our preliminary analyses show that N0 2 and CO 
did not have significant effects on mortality. 

We obtained daily data on weather conditions from the Meteorological 
Department, Ministry of Information and Communication Technology in 18 pro- 
vinces from 1999 to 2008 (Meteorological Development Bureau 2010). Daily weather 
data include four weather variables: mean, minimum, and maximum temperatures 
(°C), and relative humidity (%). 

Statistical analysis. We used two-stage Bayesian hierarchical statistical models to 
estimate province- specific and national average associations between daily air 
pollutants and mortality. In the first stage, we applied a time -stratified case-crossover 
design, commonly used to examine the effects of short-term exposures, to assess the 
associations between air pollutants and mortality in each province. The case- 
crossover design is a form of matched case-control analysis, with each case acting as 
its own control 33 . This design can successfully control for potential time -invariant 
confounding factors (e.g., smoking and obesity). When the case- crossover design is 
applied to time series data, it controls for seasonal and secular trends by matching case 
and control days in relatively small time windows (e.g., calendar month). This 
controls for season using a step function rather than a smooth spline function 34 . The 
case-crossover design fitted by conditional logistic regression (survival regression) is a 
special case of Poisson regression models 35 . This equivalence provides computational 
convenience and permits model checking for the case-crossover design using 
standard log-linear model diagnostics. 

In this study, we used a Bayesian statistical inference with Poisson regression model 
to fit a case-crossover design. It has been well acknowledged that Bayesian meth- 
odology provides enormous advantages over traditional methods for epidemiological 
analyses 36 . We used control days on the same day of the week in the same calendar 
month as the case day, to adjust for day of the week and to avoid overlap bias 37 . To 
completely control for the potential effects of ambient temperature on mortality, we 
included 15-day moving average temperature and 15-day moving average relative 
humidity with a nonparametric smooth to the models 34,38,39 . 

After establishing the basic model, we introduced the air pollutant concentrations 
into the models. We examined the associations with different lag structures, including 
single-day lag (from lag 0 to lag 5) and moving average lag (Iag0-1 to lag 0-5). For 
single-day lag models, a lag of 0 days (lag 0) corresponds to the current-day air 
pollution concentration, and a lag of 1 day (lag 1) refers to the previous day's con- 
centration; For moving average lag models, lag 0-1 corresponds to an 2-day moving 
average of air pollution concentration of the current and previous 1 days. For each 
pollutant, we fitted both single-pollutant and multiple -pollutant models to assess the 
stability of the associations. The lags producing the highest effect estimates were put 
into the multi- pollutant models. We used PM 10 as an example, in the single -pollutant 
models, PM 10 was included alone in the model; in the 2-pollutant models, PM 10 (lag 
0) and 0 3 (lag 0-5) (or S0 2 (lag 0-3)) were included jointly at the same lag. In 3- 
pollutant models, PM 10 (lag 0) and 0 3 (lag 0-5) and S0 2 (lag 0-3)were included 
jointly at the same lag. To examine whether effects of air pollutants on mortality 
differed by season, we also stratified analyses by three seasons (summer, rainy season 
and winter) in Thailand. 

In the second stage, we used Bayesian hierarchical meta-analysis to pool the 
national effect estimates of the association of air pollutants and mortality 40 . This 
method has been widely used in the area of air pollution and mortality/morbidity 41 . 
This approach provides a flexible tool to pool risk estimates while accounting for 
between -province variability (heterogeneity) and within-province statistical error of 
"true" risks. This Bayesian model produced a posterior probability distribution of 
pooled estimates, from which we reported the posterior mean and 95% posterior 
interval (PI). 
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We conducted a series of sensitivity analyses to assess the robustness of our results 
by removing temperature and relative humidity, changing time windows (14, 21 and 
28 days), introducing N0 2 and CO into the models. We also used empirical models 
(generalised linear model and generalised additive models) to check our results. 

R software (version 2.30.1) was used to perform all analyses. The "INLA" package 
was used to fit Bayesian statistical inference 42 and TLNISE package was used to fit 
Bayesian meta-analysis. The results are presented as the percentage change in daily 
mortality per 10 ug/m 3 increase of PM 10 , per 10 ppb increase of 0 3 , and per 1 ppb 
increase of S0 2 . 
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